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1 Introduction 

We summarize the discussion on the technical possi- 
bihty of the reahzation of ultracold neutron (UCN) 
source at the Japan Proton Accelerator Research 
Complex (J-PARC) for new physics beyond the stan- 
dard model [ll[2j. We discussed the UCN production 
according to the down-scattering in the superfiuid 
helium of very cold neutrons from the spallation re- 
action with the proton beam from the linear accelera- 
tor (LINAC) of the J-PARC. In this report, we report 
the summary of taskforce meetings in the period of 
Jan. -Mar. 20090. In the quickest case, J-PARC can 
provide UCNs for physics experiments after two-year 
construction period. We also expect that the UCN 
source would be the most intense source in the world 
at the startup stage and it can be reinforced by one 
order of magnitude. We recommend to consolidate 
the nation-wide and international collaboration to 
study the details of the source design, experimental 
techniques and physics experiments. 

2 Distribution of Proton Beam 

The layout of the linear accelerator (LINAC) and 
the 3 GeV Rapid Cycle Synchrotron (RCS) of the 
J-PARC is shown in Fig. [H 

The J-PARC LINAC is capable of providing pulsed 
proton beams with the repetition rate of 50 Hz, while 
it is currently providing protons to the RCS with the 
repetition rate of 25 Hz. We consider the case that 
the beam delivery to the RCS is 25 Hz without dis- 
turbing the utilization downstream of the RCS and 
the additional 25 Hz is delivered to the UCN source 
as shown in Fig. O The beam power to the UCN 
source can be adjusted by combining the following 
tunable factors. 

1. Macro Pulse Width 

The original macro pulse time width is 500/is. 
The average current value is decreased by short- 
ening this to tmacro The adjustable range is 



50/iS < tn 



< 500/is. 



(1) 



1 2009/01/28 13:00-15:00 
2009/02/04 15:30-16:30 
2009/02/12 16:00-17:00 
2009/02/19 17:00-18:00 
2009/03/11 10:30-16:00 



Chopping Ratio 

Each macro pulse is a set of intermediate pulses. 
We denote the width and the repetition rate of 
the intermediate pulse as tmicro and Tmicro- The 
standard width is (tmicro)std = 560 ns. The stan- 
dard repetition rate is chosen for the synchro- 
nization to the revolution time of protons in the 
RCS before the acceleration as 



(^micro)std 



1065ns 
815ns 



ISlMeV), 
400MeV). 



(2) 

The average beam current can be changed by 
adjusting tmicro and Tmicro- We refer to the ratio 
of the average current to that without chopping 
as the the chopping ratio. The chopping ratio is 
given as 



Trr 



"^chop 



(tmicro)std (Tmicro)std 



(3) 



The chopping ratio can be lowered by shortening 
the width and lowering the repetition rate of 
the intermediate pulse. The pulse width in the 
range of 

tm\c.r 



1 ^ 

5 (tmicro)std 



< 1 



(4) 



is already proven experimentally. On the other 
hand, the adjustable range of the repetition rate 
with the present control system is 



(Tmicro )std 2 4 



(5) 



Therefore, the proven range in the chopping ra- 
tio is 



0.05 < rchop < 1- 



(6) 



3. Switch On/ Off of Macropulses 

The average beam current can be decreased by 
switching off the beam supply from the ioninza- 
tion source. The macropulse rate fraction to the 
UCN source can be adjusted in the range of 



iVuCN < 

- 25Hz - ' 



(7) 



where we denote the macropulse rate to the 
UCN source as A^cn 




Figure 1: Layout of the LIN AC and 3GeV Rapid Cycle Synchrotron of the J-PARC. 1,2, and 3 in the figure 
are the UCN source locations discussed in the taskforce. 
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Figure 2: Pattern diagrams of beam distribution of J-PARC LIN AC. 



The peak current (/peak) is presently /peak = 
30mA with Ep = ISlMeV, which is to be upgraded 
to /peak = 50mA with Ep = 400MeV. We expect 
that the average current (/ucn) to the UCN source 
can be adjusted in the range of 



< /ucn < 



t 



^chop 



500/is 
375/iA 
625/iA 



25Hz 

iEp = 



181,400MeV) 
400MeV) 



(8) 



together with the arbitrary macropulse pattern. The 
available proton beam power to the UCN source 
W^UCN is in the following range. 



< W^UCN < 



67.9kW 
150kW 
250kW 



{Ep = ISlMeV) 
{Ep = 400MeV) 
{Ep = 400MeV) 



(9) 



In this report, we take the case that tpw = 50 and 
^UCN = 25 Hz, which leads to 

I^ucN = 20kW. (10) 



3 Proton Beam Transport 

The feasibility to transport the proton beam to the 
candidates of UCN source locations location of the 
UCN source has been discussed under following re- 
strictions. 



Curvature Radius of Beam Orbit 
The output of the LINAC is ions. The 
strength of the magnetic field to bend the H" 
beam is restricted to be less than 0.5T to avoid 
the strip-off of the electrons. 



An alternative solution is to convert the charge 
state from H~ to 11+ by installing a charge strip- 
ping foil, which suppresses the orbit curvature 
radius by a factor of 2. However, we assume 
to use the H~ beam to avoid additional efforts 
to maintain the stripping foil and the radiation 
shields. 



• Beam Size 

The lifetime of the proton incident window is in- 
versely proportional to the proton current den- 
sity as Eq. [m The proton current density needs 
to be appropriately controlled for a reasonable 
operation lifetime since the proton incident win- 
dow of the neutron production target experi- 
ences the most severe radiation damage. We 
require that the beam size is expanded to be 
about 10 ~ 25mm (r.m.s.) on the proton inci- 
dent window. 



3.1 UCN Location Candidate 1: Ex- 
tension of the LINAC down- 
stream to the south 

Here we discuss the case of the UCN source is lo- 
cated at the position 1 in the Fig. [TJ The proton 
beam is branched as shown in Fig. [3l The proton 
beam is kicked off to the west and transported off 
the original LINAC line so that it does not disturb 
the "Originally planned SCL R&D line" of Fig. H 
which is temporarily reserved for the future exten- 
sion to the 2nd neutron target station. We assume 
that the UCN source would be installed in the ex- 
tended building on the same level as the floor level 
of the LINAC building. 

The proton beam is transported for about 60 m 
with the angle of 4o off the LINAC axis. The proton 
beam optics is shown in Fig. [H The beam size at the 
target position is Xrms — 20mm and yrms — 10mm. 
The window lifetime is estimated to be about 27 
years according to Eq. [TH with the operation time 
of 5000 hours per year even for Ep = ISlMeV case, 
which is long enough in terms of the maintenance 
efforts. The relation between the UCN source and 
the "Originally planned SCL R&D line" is shown in 
Fig. [5l The width for the LINAC extension is re- 
stricted in about 10 m long region. The size of the 
neutron source is described later (c.f. section [4^ . 

3.2 UCN Location candidate 2: Ex- 
tension of the LINAC to south- 
west 

It is possible to transport the proton beam further 
off the LINAC axis to the location 2 in Fig. [1] by 
additionally bending to the west as shown in Fig. [6l 
However, this candidate has not been discussed in 
details. 

3.3 Location Candidate 3: U-turn 
Layout 

The location 3 in the Fig. [T] is also possible by bend- 
ing the proton beam by 180°. 

In this case, the first 90° bending has the same 
curvature as the L3BT transport line in the existing 
tunnel. The beam is branched from the ADS line and 
bent by another 90°. The beam defocus is initiated 
after the beam branch from the ADS line so that it 
does not affect the beam optics to the ADS line. 

The quadrupole magnets are not installed in the 
second 90° bending for minimizing the curvature ra- 
dius since the achromat beam transport is not neces- 
sary for the UCN line. The curvature radius of 8.3 m 
can be realized by employing six magnets of 1.67 m 
long, 15° bending angle with 6.37 m curvature radius 
placed with the spacing of 0.6 m. 

However, the drift space for expansion of the beam 



Ep 


(^UCN)max 


proton kinetic energy 


beam current 




for 20 kW operation 


181 MeV 


110/iA 


400 MeV 


50/iA 



Table 1: Beam current for 20kW beam power. 



size can be decreased by installing two or three 
quadrupole magnets during the 90° bending. The 
length of the drift space is estimated as 20 m down- 
stream of the 90° bending. 

The corresponding UCN source location is shown 
in Fig. [71 0. In this design, the beam size on the 
neutron production target is 2xrms ~ 15mm and 
2^rms ~ 15mm. The beam size corresponds to the 
window lifetime of about 25 years, which is reason- 
ably long. 

4 Ultracold Neutron Source 

We assume the proton beam with the kinetic energy 
of 181 MeV or 400 MeV is incident to the neutron 
production target with the beam power of 20 kW. 
Neutrons are generated by the spallation reaction 
and moderated into the cold region. The cold neu- 
trons are converted to ultracold neutrons according 
to the inelastic scattering in superfluid helium. We 
also assume that ultracold neutrons are accumulated 
in the superfluid helium. 

4.1 Production Target 

The range of a 400 MeV proton in Pb is about 13 cm 
as shown in Fig. [9] (simulation code is PHITS [3j). 
We, hereafter, adopt the use of a 20 cm thick Ta- 
coated W target, which has been established as the 
spallation neutron target material. 

The proton incident window of the neutron pro- 
duction target experiences the most severe radiation 
damage. The lifetime of the window is estimated by 
evaluating the displacement per atom (dpa) H. 

We estimate the case of aluminum alloy A5083 
which is used in the J-PARC MLF neutron source. 
The dpa rate is estimated as dpa = 0.44 y~^ for the 
peak current density of 4.3/iA cm~^ with the 5000 
hour operation in a year [4J. If we assume the ac- 
ceptable dpa is (dpa)max = 10, the lifetime is given 

^ If we use H+ beam, the curvature radius can be half since 
the magnetic field can be twice. However, the drift space 
required for the beam defocus doesn't change and only the 
curvature of the second 90° bending can be shortened. We 
adopt H~ beam transport considering both of the advantage 
and disadvantage. 

^ Although an additional effect of brittleness caused by He 
and H from the nuclear reaction induced by 3 GeV proton 
beam, the dpa evaluation is considered appropriate for 400 
MeV protons. 
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Figure 3: UCN source location candidate 1. 
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Figure 4: Proton beam optics to the UCN source location candidate 1. 
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Figure 5: Relation between the UCN source and the "Originahy planned SCL R&D line". 




Figure 6: UCN source location candidate 2. 
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Figure 7: UCN source location candidate 3. 
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Figure 8: Proton beam optics to UCN source location candidate 3. 2xrms and 2?/rms are shown. 
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Figure 9: Distribution of 400 MeV protons in the 
Pb target. 

as4r=3GeV>22.7y0. 

One of the difference between our case and J- 
PARC MLF neutron source is the proton energy. 
The Fig. 2 in the reference [4J shows the displacement 
cross section at = 181 MeV, 400 MeV and 3 
GeV are crudely related as 



cTd{Ep = ISllMeV]) 



cTdiEp = 3 [GeV]) X 1.3. 

(11) 



Another difference is the beam distribution. While 
the proton beam position is scanned to paint the win- 
dow area so that the average beam density becomes 
as uniform as possible for the J-PARC MLF neu- 
tron source, the beam distribution can be approxi- 
mated by a Gaussian distribution for our case. We 
assume the root-mean-square along x-axis (horizon- 
tal) and that along y-ax.is (vertical) as Xrms and ^rms- 
We put the average current density in the region of 

-2Xrms < X < 2Xrms and -2?/rms < V < 2y 

rms as Jo- 

The current density at the beam center is given as 
ji'c 2. 8 jo and can be written as 



jc - 2.^ 



4 



IGx^jxis^ri 



(12) 



where /ave is the average current. The average cur- 
rent is /ave = 1.1 X 10~^ A and /ave = 5 X 10~^ A for 
Ep = 181MeV and Ep = 400MeV if we assume the 
average beam power of VFucn = 20 kW. 

Therefore, we can estimate the window lifetime as 



,dpa=10 
^£;„ = 181MeV 



dpa=io cTdjEp 3 [GeV]) 
^E,=3GeV ^^^^^ = 181 [MeV]) 

4.3/iAcm~^ 



2.8/ave/ (ISXrnis^rms) 



^ There is another estimation with the peak current of 
10/iA cm~^ and (dpa)ma 
11.3 y. 



5, which leads to ^^^t=3^eV — 



10mm 10mm 

^rms ^rms 

25mm 25mm 



X 11 year 
X 68year 



(13) 



dpa=10 



=400MeV 



dpa^io cidjEp = 3[GeY]) 
=3^^^arf(£;p = 400[MeV]) 
4.3/iAcm~^ 



2.8/ave/(16cra,Cr^) 
^rms Vrms 



10mm 10mm 

^rms ^rms 

25mm 25mm 



X 24year 
X 150year (14) 



4.2 Neutron Moderator 

We consider the conversion of cold neutrons into ul- 
tracold region (/^ucN = 252 neV) out of thermal 
equilibrium according to the inelastic scatteing. An 
intense radiation field of neutrons as cold as pos- 
sible is necessary to increase the ultracold neutron 
intensity. A 20 kW cold source provides the neu- 
tron fiux of about $0 — 2 X 10^^cm~^s~^, which 
is crudely estimated assuming the 400 MeV protons 
are incident to a Pb target with the average current 
of 50/iA and moderated in D2 with the refiector of 
Be and graphite as shown in Fig [141 THe neutron 
temperature is about 30K. 

The ultracold neutron density accumulated in the 
converter can be given as 



PUCN = Pr{l - e-'l^) 



(15) 



where P is the conversion rate, r storage time. The 
values of P and r are listed in table [2] [6j . The time 
evolution of Eq. [15] using the values in table [2] is 
shown in Fig.[TOl The UCN density strongly depends 
on the volume ratio and geometrical arrangement of 
converter volume and storage volume. However, in 
this summary, we assume that the converter volume 
is the same as the storage volume. Therefore, the 
value is correct as it is, and remaining corrections are 
the UCN loss on the storage wall, due to upscattering 
and absorption in the converter. 

The reversal process of the super-cold neutron gen- 
eration happens in the superfiuidity helium, too. It 
happens by short odds by calling scattering (upscat- 
teirng) on, and the phonon density high (By the tem- 
perature of the superfiuiditive helium high) as for 
the reversal process. As for the constant, when a 
super-cold neutron is lost by the above one scatter- 
ing, what shown as a function of the temperature 
of the superfiuiditive helium is Fig. [Til Because the 
time constant of the loss by the above one scatter- 
ing is a neutron longevity level if the temperature 
of the superfiuiditive helium is kept about less than 
0.8K, the loss by the above one scattering can be 
disregarded 0. 



^ Because the neutron absorption by ^He that exists nat- 
urally by about Ippm cannot be disregarded when actually 



converter 


He-II 


Solid ortlio-D2 


a-02 


excitation 


phonon 


phonon 


magnon 


converter temperature 


0.7 K 


5 K 


2 K 


appropriate neutron temperature 




on 


1 o ly 
iz K 


conversion rate 








Icm-Sg-i <|)q 

@(cold neutron temperature of 30K) 


1 Qon 






storage time r 


886 s 


0.146 s 


0.750 s 


@(wall loss and up-scatteirng ignored) 








maximum UCN density in the converter 


1,600,000 


3,800 


21,000 









Table 2: Characteristics of UCN converter materials [6 . The UCN density strongly depends on the volume 
ratio and geometrical arrangement of converter volume and storage volume. However, in this summary, we 
assume that the converter volume is the same as the storage volume. Therefore, the value is correct as it is, 
and remaining corrections are the UCN loss on the storage wall, due to upscattering and absorption in the 
converter. 
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Figure 10: Time evolution of the density of ultracold 
neutrons accumulated in the converter assuming that 
the temperature of cold neutrons is 30 K (Boltzmann 
distribution) and the flux is <l>o = 2 x 10^^ cm~^ s~^. 
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Figure 11: Temperature dependence of the time con- 
stant of upscattering in the superfluid helium (He- 
H). 



4.3 Heat Load 

The incident beam would be tailored by a tapered 
collimator so that the central 90% beam power 
reaches to the neutron production target. The neu- 
tron production target is the Ta-coated W. The tar- 
get thickness is 20 cm and the cross section is almost 
as large as the beam cross section. About 10% of the 
beam power is consumed in the spallation reaction. 
The distribution of the heat load is estimated follow- 
ing the design of the J-PARC MLF neutron source 
as shown in table [3] and Fig. [121 

The heat load in the neutron production target 
can be removed by the water cooling and the solid 
target can be used. 

The heat load in the moderator is sufficiently small 



accumulating in the superfluiditive helium, it is necessary to 
use super-high purity (concentrate ^He). 





heat load [W] 


fraction 


collimator 


2000 


0.1 


neutron production target 


9000 


0.45 


moderator 


72 


0.0036 


UCN converter 


10 




reflector and shield 


7200 


0.36 



Table 3: Distribution of the heat load. Missing heat 
is consumed in the spallation reaction. 



proton beam 




Figure 12: Distribution of the heat load. 

to be removed by the evaporation of liquid helium. 
Pressurized circulation should be considered beyond 
300 W heat load when the beam power is increased. 

The heat load in the UCN converter is about 10 W. 
The required cooling power is estimated in terms of 
the consumption of liquid helium at several converter 
temperatures. 

• (1) 5K 

The helium consumption is 15 L/h and 360 
L/day, which is calculated from the latent heat 
of liquid helium (20 kJ kg-^). A 1000 L vessel 
should be supplied every two days. 

• (2) 2K 

The helium consumption is 25 L/h and 600 
L/day, which is estimated based on the experi- 
ence of the Low Temperature Engineering Cen- 
ter of KEK: 15 W cooling power at 2 K. A 1000 
L vessel should be supplied every 1.3 days. 

• (3) 0.8K 

The helium consumption is 50-100 L/h and 
1200-2400 L/day, which is a crude estimation 
since we have no experience. One to three 1000 
L vessel should be supplied every day. 

In the present discussion, the UCN converter is as- 
sumed to be kept at 0.8 K under the 15 W heat load. 

A new liquefier or refrigerator is recommended to 
satisfy the corresponding consumption of liquid he- 
lium |j. A new liquefier requires the increase of qual- 
ified manpowers for the safety and operation. If we 

^ The practical way to employ a dedicated liquefier is to 



shield material 


attenuation length 
(3 GeV proton) 


attenulation length 
(400 MeV proton) 


iron 
concrete 


0.244 m 
0.650 m 


0.177 m 
0.409 m 



Table 4: Attenuation length of iron and concrete |5]. 



directly connect the refrigerator and the UCN con- 
verter, we may decrease the amount of qualified man- 
power for safety and operation spending more time 
and efforts to get the admission by the Government. 

4.4 Radiation Shield 

We require that the radiation dose rate on the shield 
surface is 1/iSv h~^ or less, consistently with the 
safety r^ulation of the J-PARC MLF experimen- 
tal hallQ This safety regulation was determined for 
the safety and also for the reasonably low radiation 
background for physics measurements. 

We determined the dimension of radiation shield 
by scaling the J-PARC MLF neutron source by the 
beam power. 




Figure 13: Radiation shield structure used in the 
design. 

In the case of the UCN location candidate 1, the 
UCN source is neighboring to the tunnel for the 
future extension of the linear accelerator, in which 
there is no access during the beam operation. If we 
assume that the radiation dose rate of ImSv h~^ is 
acceptable during the operation in the accelerator 
tunnel, a 1.5 m thick iron wall is sufficient for the 
beam power of 20 kW. 

Under these assumptions, the minimum radiation 
shield is 117m^ (916 tons) of iron and 147m^ (324 
tons) of concrete as shown in Fig. [HI The rectangu- 
lar hexahedron tangential to the minimum column 
is 150m^ (1150 ton) of iron and 169m^ (373 ton) of 
concrete. 

separate it from the UCN converter, which consumes the liq- 
uid helium, without connecting the liquefier and the UCN 
converter. Otherwise, the national law on the safety of high 
pressure gas requires a special admission by the Minister of 
Industry of the Japanese Government. 

^ The design goal of the J-PARC Controlled Area Category 
1 is 12.5/iSv h"-*^. This radiation dose rate is also considered 
and the result is shown in Table [5] 





MLF 


UCN 


UCN 


UCN 


UCN 


Ep [MeV] 


3000 


400 


400 


400 


400 


Power [kW] 


1000 


20 


20 


200 


200 


surface dose rate 


1 


1 


12.5 


1 


12.5 


@[/iSv h-i] 












core radius i^xMRA [m] () 


0.687 


0.69 


0.69 


0.69 


0.69 


iron thickness I^iron [m] 


4.113 


2.6 


2.4 


2.8 


2.6 


concrete thickness I^conc [^] 


3.200 


1.2 


0.7 


1.7 


1.2 


outer diameter [m] 


16 


9.1 


6.2 


10.4 


8.9 



Table 5: Estimation of required size of the radiation shield. 
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Figure 14: Configuration and outer dimension of the 20 kW UCN source with the neutron production target, 
moderator, refiector, converter and radiation shield with the surface dose rate of 1/iSv/h. In this example, the 
radiation shield on the left rear is decreased assuming that there is no access to that side during the operation. 



5 Construct Cost 



6 Near Fugure 



The dimension of the building to instah the UCN 
source is shown in Fig. [151 The floor level of the 




Figure 15: Building for the UCN source and the ex- 
periment for the case of UCN source location candi- 
date 1. The tunnel for the future extension of the lin- 
ear accelerator is on the left rear of the UCN source. 



building is 13m below the ground level which is the 
accelerator building floor level. The building dimen- 
sion is d' = 5m a = 30m, b = 27m, h = 8m, and 
d = 13m as shown in Fig. [161 The construction cost 
of the building is estimated as about 950 million yen 
in rough figures. 




Details of the design and the physics program using 
the UCN source will be discussed in the fiscal year 
of 2009, which consequently expand the researchers' 
community. Measurement techniques will be stud- 
ied and developed using the cold beam line for the 
neutron optics at the BL05 of J-PARC MLF. A new 
research collaboration framework will be promoted 
through above activities at nation-wide level and also 
at the international level. The collaboration frame- 
work will pay efforts to obtain sufficient financial sup- 
port. 
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Figure 16: Definition of the outer dimension of the 
building. 



The rough cost estimate for the UCN source con- 
struction is listed in Table [6l The construction 
cost can be suppressed by the re-use of proton beam 
transport and radiation shield and by careful re- 
design of the building. The period of the construc- 
tion is estimated as 18 months without any contin- 
gency. 





(million yen) 




building construction 


950 




proton beam branch 


pulse magnet 


50 




and transport 


DC magnets 


25 






quadrupoles and power supplies 


15 






control 


25 






linear part quadrupoles 


25 


operated by 




linear part quadrupole power supplies 


20 


present LINAC team 




steering magnet 


10 






steering magnet power supply 


10 






vacuum 


25 






beamline radiation shield 


15 




UCN source 


production target 


45 


3 personnel 




moderator 


100 


(1 personnel all the time) 




converter 


120 






reflector and radiation shield 


150 




heium liquefier 


(120L/h) 


400 


1 qualified personnel 








during operation 


neutron transport 




1 technician 


physics experiments 




not involved in this report 


total 


1,985 





Table 6: Crude cost estimation of the UCN source construction at the J-PARC. 
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